The structure of Green functions in quantum 
field theory with a general state 

Christian Brouder 

Abstract. In quantum field theory, the Green function is usually calculated as 
the expectation value of the time-ordered product of fields over the vacuum. In 
some cases, especially in degenerate systems, expectation values over general 
states are required. The corresponding Green functions are essentially more 
complex than in the vacuum, because they cannot be written in terms of 
standard Feynman diagrams. Here, a method is proposed to determine the 
structure of these Green functions and to derive nonperturbative equations 
for them. The main idea is to transform the cumulants describing correlations 
into interaction terms. 
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1. Introduction 

High-energy physics uses quantum field theory mainly to describe scattering ex- 
periments through the S-matrix. In solid-state or molecular physics, we are rather 
interested in the value of physical observables, such as the charge and current 
densities inside the sample or the response to an external perturbation. At the 
quantum field theory (QFT) level, these quantities are calculated as expectation 
values of Heisenberg operators. For example, the current density for a system in a 
state 1$) is ($|J(a;)|$), where |$) and J(a::) are written in the Heisenberg picture. 

The first QFT calculation of Heisenberg operators was made by Dyson in 
two difficult papers [TJ [5] that were completely ignored. At about the same time, 
Gell-Mann and Low discovered that, when the initial state of the system is non- 
degenerate, the expectation value of a Heisenberg operators can be obtained by a 
relatively simple formula [31 . The Gell-Mann and Low formula has been immensely 
successful and is a key element of the many-body theory of condensed matter [4 , 5 . 
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Its main advantage over the formalism developed by Dyson is that all the standard 
tools of QFT can be used without change. 

However, it was soon realized that the assumption of a nondegenerate initial 
state is not always valid. As a matter of fact, the problem of what happens when the 
initial state is not trivial is so natural that it was discussed in many fields of physics: 
statistical physics [6 , many-body physics [7] , solid-state physics , atomic physics 
[5j, quantum field theory and nuclear physics [lOllllj . As a consequence, the theory 
developed to solve this problem received several names such as nonequilibrium 
quantum field theory (or quantum statistical mechanics) with initial correlations 
(or with cumulants, or for open shells, or for degenerate systems). It is also called 
the closed-time path or the (Schwinger-)Keldysh approach for an arbitrary initial 
density matrix. 

It should be stressed that the problem of the quantum field theory of a 
degenerate system is not only of academic interest. For instance, many strongly- 
correlated systems contain open-shell transition metal ions which are degenerate 
by symmetry. This degeneracy makes the system very sensitive to external per- 
turbation and, therefore, quite useful for the design of functional materials. 

The elaboration of a QFT for degenerate systems took a long time. It started 
with Symanzik 12J and Schwinger [13^ and made slow progress because the com- 
binatorial complexity is much higher than with standard QFT. To illustrate this 
crucial point, it is important to consider an example. According to Wick's the- 
orem, the time-ordered product of free fields can be written in terms of normal 
order products: 

T(p{xi) . . .ip{x4) = -.ifixi) . . .ip{x4): + '^:ip{xi)ip{xj):Go{xk,xi) 

ijkl 

+ y^/.tf{xk)(p{xi): Go{ j)Go{xk,xi), 

ijkl ijkl 

where the quadruplet of indices (i, j, fc, I) runs over (1, 2, 3, 4), (1, 3, 2, 4) and (1, 4, 2, 3). 
The expectation value of this expression over the vacuum gives the familiar result 
Syfci G'o(a;i, a;j)G'o(xfe, xj). However, when the initial state is not the vacuum 
(as in solid-state physics), we obtain 

{ip\T(p{xi) . . . Lp{x4)\il;) = {ipl-ifiixi) . . . ifi{x4):\ip) + ^ p2{xi, Xj)Go{xk, xi) 

ijkl 

+ X! P^i^k, xi)Go{xi, Xj) + E Go{xi,Xj)Go{xk,xi), 

ijkl ijkl 

where p2{x,y) = {tlj\:i^{x)(p{y):\'ip) . If we assume, for notational convenience, that 
the expectation value of the normal product of an odd number of field operators 
is zero, the fourth cumulant P4{xi, . . . , x^) is defined by the equation 

{ipl-ip^Xi) . . .ip{x4,):\ljj) = pi{xi, . . . ,X4) +^p2{Xk,Xl)p2{Xi,Xj). 

ijkl 
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If we put g = Gq + p2, the free four-point Green function becomes 

{tp\Tip{xi)...Lp{x4)\ip) = p4,{xi,...,X4) + ^g{xi,Xj)g{xk,xi). 

ijkl 

When /94 — 0, the expression is the same as over the vacuum, except for the fact 
that the free Feynman propagator Go is replaced by g. When this substitution is 
vahd, standard QFT can be apphed without major change and the structure of 
the interacting Green functions is not modified. For fermionic systems described 
by a quadratic Hamiltonian Hq-, this happens when the ground state is nondegen- 
erate, so that \ip) is a Slater determinant. When pi ^ 0, the expression becomes 
essentially different because the cumulant p4 appears as a sort of free Feynman 
propagator with four legs. In general, the expectation value of a time-ordered 
product of n free fields involves pk with k < n. 

In other words, the perturbative expansion of the Green functions can no 
longer be written as a sum of standard Feynman diagrams. Generalized Feynman 
diagrams have to be used, involving free Feynman propagators with any number 
of legs [ilZlITl]. 

Because of this additional complexity, the structure of the Green functions 
for degenerate systems is almost completely unknown. The only result available is 
the equivalent of the Dyson equation for the one-body Green function G{x,y) [7] 

G = (l-A)-i(Go + G)(l-S)-i(l + SG), 

where A, _B, C and S are sums of one-particle irreducible diagrams. When the 
initial state is nondegenerate, A = B = G = and the Dyson equation G = 
Go + Go EG is recovered. 

In the present paper, a formal method is presented to determine the structure 
of Green functions for degenerate systems. The main idea is to use external sources 
that transform the additional propagators /?„ into interaction terms. This brings 
the problem back into the standard QFT scheme, where many structural results 
are available. 

2. Expectation value of Heisenberg operators 

Let us consider a physical observable for instance the charge density or the 

local magnetic field. In the Heisenberg picture, this observable is represented by 
the operator Ah {t) and the value of its observable when the system is in the state 
<i>_ff} is given by the expectation value {A{t)) = 

Going over to the interaction picture, we write the Hamiltonian of the system 
as the sum of a free and an interaction parts: Ii{t) — Hq + Hj{t), we define the 
evolution operator U{t,t') = T(exp(— i J^, Hi{t)dt)) and we assume that the state 
\^h) can be obtained as the adiabatic evolution of an eigenstate |$o} of Hq. The 
expectation value of A becomes 

{A{t)) = {<!>o\U{~oo,t)A{t)U{t,^^)\<i>o), 
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where A{t) on the right hand side is the operator representing the observable 
in the interaction picture. The identity 1 = U{t, oo)U{oo,t) and the definition 
S = U{oo, — oo) enable us to derive the basic expression for the expectation value 
of an observable in the interaction picture: 

{A{t)) = {^o\S^T{A{t)S)\'^>„). (2.1) 

When |$o} is nondegenerate, this expression can be further simplified into the 
Gell-Mann and Low formula 

If the system is in a mixed state, as is the case for a degenerate system by 
Liiders' principle, the expectation value becomes 

{A{t)) = 5]p„(<i>„|5tr(^(i)5)|$„)^ 

n 

where Pn is the probability to find the system in the eigenstate \^n)- It will be 

convenient to use more general mixed states J2mn^mn\^m) {^n\, where 

a density matrix (i.e. a nonnegative Hermitian matrix with unit trace). Such a 

mixed state corresponds to a linear form ui defined by its value over an operator 

O: 

m/ • 

mn 

Then, the expectation value of A{t) becomes 

{A{t)) = u;{S^TiAit)S)). (2.2) 



3. QFT with a general state 

In all practical cases, the operator representing the observable A{t) in the interac- 
tion picture is a polynomial in ip and its derivatives. Its expectation value (|2.2p can 
be expressed in terms of Green functions that are conveniently calculated by a for- 
mal trick due to Symanzik [12J and Schwinger [13], and reinterpreted by Keldysh 

The first step is to define an S-matrix in the presence of an external current 
j as S{j) = T(e-*/^""(*)^*+'/j(^)'^(^)'^="), where i?'"* in the interaction Hamil- 
tonian in the interaction picture. The interaction Hamiltonian is then written in 
terms of a Hamiltonian density V{x), so that J H"^^{t)dt = J V{x)dx and the 
generating function of the interacting Green functions is defined by Z{j+,j-) = 
uj(^S^ {j-)S{j+)) . The interacting Green functions can then be obtained as func- 
tional derivatives of Z with respect to the external currents j+ and j_ . For example 

(T(v.(x)y>(y)))=- /.'f^+'.^;\ , and {^(a^My)) = fl^i±ll-l- 

As in standard QFT, the connected Green functions are generated by log Z. 
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In the functional method [HIIIJ, the generating function Z of the interacting 
system is written as Z = e~*^Zo, where D is the interaction in terms of functional 
derivatives 

and where Zoij+,j-) - Lu{Sl{j_)So{j+)) , with So{j) = r(e'/j(-)'^(-)d-). Note 
that Zo(j+j J-) is the generating function of the free Green functions. 
A straightforward calculation fTTI leads to 

= e^-^/^ •^■'''^^'^o^'^'^'j'^^^'^'*^e'''^^+~^^-', 

where j — (j+ii-) the source vector, 

^ _ MO|T(0(a:)0(2/))|O) -(O|</.(y)0(a:)|O) \ 

- 1, -(O|0(a;)0(y)|O) (O|T(0(x)0(y))|O) j' ^^'^^ 

is a free Green function (with T the anti-time ordering operator) and 

defines the generating function p'{j) of the cumulants of the initial state w. 

The free Green function G[, describes the dynamics generated by the free 
Hamiltonian Hq. It can also be written in terms of advanced and retarded Green 
functions [13] . 

The idea of describing a state by its cumulants was introduced in QFT by 
Fujita ^ and Hall ^2,. It was recently rediscovered in nuclear physics [TUl [TT| and 
in quantum chemistry [18] . 

The next step is to modify the definition of the free Green function. The 
cumulant function is Taylor expanded 

P'(^) = dXl...dXnPn{Xl,...,Xn)j{Xl)...j{Xn)■ 
The expansion starts at n = 2 because — 1 and the linear term can be 
removed by shifting the field ip. The bilinear term p2{x, y) is included into the free 
Green function by defining 

GQ{x,y) = G'f:,{x,y) ^ p2{x,y)(^ 
and the corresponding cumulant function becomes 

P{j) = P'i.i) - {'^/^) J dxdyjix)p2{x,y)j{y) 

= dXi...dXnPnixi,...,X„)j{xi)...j{Xn)- 

n— 3 
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Remark 3.1. There are several good reasons to use Go and p instead of G'q and p': 
(i) This modification is exactly what is done in solid-state physics when the free 
Green function includes a sum over occupied states 19|; (ii) At a fundamental level, 
Go and p have a more intrinsic meaning than G'q and p' because they do not depend 
on the state |0) chosen as the vacuum; (iii) An important theorem of quantum field 
theory f20' states that, under quite general conditions, p„(a;i, . . . smooth 
function of its arguments when n > 2, so that Go gathers all possible singular 
terms (a related result was obtained by Tikhodeev [21]); (iv) A state for which 
p{j) = is called a quasi-free state 22 , quasi-free states are very convenient in 
practice because the rules of standard QFT can be used without basic changes. 
Thus, the additional complications arise precisely when p (and not p') is not zero. 



4. Nonperturbative equations 

To size up the combinatorial complexity due to the presence of a non-zero p, we 
present the diagrammatic expansion of the one-body Green function G{x, y) for 
the Lp^ theory to second order in perturbation theory. For this illustrative purpose, 
it will be enough to say that the cumulant pn{xi, . . . ,Xn) is pictured as a white 
vertex with n edges attached to it, the other vertex of the edge is associated with 
one of the points xi, . . . , Xn- For example, pi{xi, . . . , x^) is represented by the 
diagram 



In this diagram, the white dot does not stand for a spacetime point, it just 
indicates that the points xi to x^ are arguments of a common cumulant. If we 
restrict the calculation to the case when p„ = if n is odd, we obtain the following 
expansion 
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In standard QFT, only the first and last diagrams of the right hand side are 
present. In the general case when all p„ ^ 0, the number of diagrams is still much 
larger. 

4.1. Generalized Dyson equation 

As mentionned in the introduction, the only known result concerning the structure 
of Green functions with a general state was derived by Hall for the one-body Green 
function G(x, y) [7] 

G = (1 -^)-i(Go + C)(l -S)-i(l + SG). 

In diagrammatic terms the quantities A, i?, C and E are sums of one-particle 
irreducible diagrams. If we take our example of the Green function of Lp^ theory 
up to second order, we find 



S = 1-1- <=oo -I- 1 -I- o + • • • 

In standard QFT, we have A = B = C = and the diagrammatic represen- 
tation of S contains much less terms. However, the difference with standard QFT 
is not only limited to the number of diagrams. The definition (|3.2p of the cumu- 
lant function, and the fact that the free field (/? is a solution of the Klein-Gordon 
equation imply that p„ is a solution of the Klein-Gordon equation in each of its 
variables. Thus, A{x,y), B{x,y) and C{x,y) are solutions of the Klein-Gordon 
equation for x and y. As a consequence, applying the Klein-Gordon operator to 
the Green function gives us (□-|-m^)G — (1 — i3)^^(l-|-SG). In other words, apply- 
ing the Klein-Gordon operator kills a large number of terms of G. This is in stark 
contrast with standard QFT, where (□-|-to^)G = 1-1- EG and amputating a Green 
function does not modify its structure. This important difference makes some tools 
of standard QFT (e.g. amputated diagrams or Legendre transformation) invalid 
in the presence of a general state. 

All those difficulties explain the scarcity of results available in non-perturbative 
QFT with a general state. Apart from Hall's work [7J, the only non-perturbative 
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results are Tikhodeev's cancellation theorems [23l [24] and the equation of motion 
for the Green functions [25 . 

In the next section, we present a simple trick to derive the structure of Green 
functions with a general state. 

4.2. Quadrupling the sources 

We first determine the main formal difference between standard QFT and QFT 
with a general state. In both cases, the generating function of the Green functions 
can be written Z = e~*^Zo, where D describes the interaction and the initial 
state. In the presence of a general state, the interaction D is simple but is made 
non standard by the cumulant factor if . The idea of the solution is to transfer 
the cumulant function p from Zq to D, because powerful functional methods were 
developed to deal with general interactions D. These methods were first proposed 
by Dominicis and Englert [IB] and greatly expanded by the Soviet school [571 dHl 

This transfer from the initial state to the interaction can be done easily by 
introducing two additional external sources and fc_ and using the identity 



gPO+-j-) — g''^ '^'ik_)pi f{j+{x)k+{x)-j-(x)k-{x))dx\ 



The term involving p can now be transferred from Zq to D by defining the new 
generating function 

where the modified interaction is 

and the modified free generating function is 

Zoij±,k±) - e-l/2j3ix)Go{x,y)3iv)dxdy^ 

with J = (j+ ,j^,k^,k^). The modified free Green function Go is now a 4x4 matrix 
that can be written as a 2x2 matrix of 2x2 matrices 



D = V 



Gn = 



Go -il 

a 



In contrast to the standard case, the free Green function Go is invertible 

il 



'^o ^ ^ il Go 

and it is again possible to use amputated diagrams and Legendre transformations. 
The free generating function Zq is the exponential of a function that is bilinear in 
the sources, and all the standard structural tools of QFT are available again. We 
illustrate this by recovering Hall's analogue of the Dyson equation. 
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4.3. An algebrEiic proof of Hall's equation 

The free generating function Zq has a standard form and the Dyson equation holds 
again: G = Go + GqSG, where G is the 4x4 one-body Green function obtained 
from the generating function Z and E is the corresponding self-energy. Each 4x4 
matrix is written as a 2x2 matrix of 2x2 matrices. For example 

G — ( \ 
\^ G21 G22 / 

We want to determine the structure of the 2x2 Green function G, which is equal 
to Gil when ~ = 0. 

The upper-left component of the Dyson equation for G is 

Gil = Go + (GoSii-zS2i)Gii + (GoSi2-zS22)G2i. (4.1) 

The lower-left component gives us G21 = — i(l -I- iSi2)~^{l + SuGn). If we in- 
troduce this expression for G21 into equation (|4.ip . rearrange a bit and use the 
operator identity 1 + 0(1 — O)^^ = (1 — O)^^ , we obtain 

(H-iS2i)Gii = (Go-S22)(l-HzSi2)"i(l-HSiiGii). 

Hall's equation is recovered by identifying A — — iS2i, B — — iSi2 and G = — S22, 
where the right hand side is taken at k+ = k- = 0. Note that Hall's equation 
is now obtained after a few lines of algebra instead of a subtle analysis of the 
graphical structure of the diagrams. 

With the same approach, all the nonperturbative methods used in solid- 
state physics, such as the GW approximation 36J and the Bethe-Salpeter equation 
[371 [38] . can be transposed to the case of a general initial state. This will be 
presented in a forthcoming publication. 

5. Determination of the ground state 

QFT with a general state was studied because the initial eigenstate of a quantum 
system is sometimes degenerate. However, it remains to determine which density 
matrix Umn of the free Hamiltonian leads to the ground state of the interacting 
system. 

A solution to this problem was inspired by quantum chemistry methods [391 . 
A number of eigenstates of Hq are chosen, for example the complete list of 
degenerate eigenstates corresponding to a given energy. These eigenstates span the 
so-called model space and the ground state of the interacting system is assumed to 
belong to the adiabatic evolution of the model space. This model space generates, 
for each density matrix, a linear form lo as described in equation (12. 2p . The problem 
boils down to the determination of the density matrix ujmn that minimizes the 
energy of the interacting system. 

This minimization leads to an effective Hamiltonian and the proper density 
matrix is obtained by diagonalizing the effective Hamiltonian. This type of method 
is typical of atomic and molecular physics [40,. However, the effective Hamiltonian 
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can now be determined by powerful non-perturbative Green function methods. 
Therefore, the present approach leads to a sort of unification of quantum chemistry 
and QFT: it contains standard QFT when the dimension of the model space is one, 
it contains standard quantum chemistry (more precisely many-body perturbation 
theory) when the Green functions are expanded perturbatively. 

Therefore, the present approach might help developing some new nonpertur- 
bative methods in quantum chemistry. On the other hand, quantum chemistry 
has accumulated an impressive body of results. The physics Nobel-prize winner 
Kenneth Wilson stated that [41j "Ab initio quantum chemistry is an emerging 
computational area that is fifty years ahead of lattice gauge theory." Therefore, 
the experience gained in quantum chemistry can be used to solve some of the 
remaining problems of the present approach, such as the removal of the secular 
terms to all order. 



6. Conclusion 

The present paper sketched a new method to determine the Green functions of 
quantum field theory with a general state. The main idea is to transform the 
cumulant function describing the intial state into an interaction term. As a conse- 
quence, the cumulants become dressed by the interaction, providing a much better 
description of the correlation in the system. 

An alternative method would be to work at the operator level, as was done 
recently by Diitsch and Fredenhagen [42] , and to take the expectation value at the 
end of the calculation. This would have the obvious advantage of dealing with a 
fully rigorous theory. However, we would loose the non-perturbative aspects of the 
present approach. 

Although this approach seems promissing, much remains to be done before 
it can be applied to realistic systems: (i) our description is purely formal; (ii) the 
degenerate initial eigenstates lead to secular terms that must be removed |T4]; 
(iii) renormalization must be included, although this will probably not be very 
different from the standard case, because all the singularities of the free system 
are restricted to Gq. 

Interesting connections can be made with other problems. For example, the 
cancellation theorem [23] seems to be interpretable as a consequence of the uni- 
tarity of the S-matrix. It would extend Veltman's largest time equation [43] to the 
case of spacetime points with equal time. Another exciting track would be a con- 
nection with noncommutative geometry. Keldysh[T5] noticed that the doubling of 
sources could be replaced by a doubling of spacetime points. In other words, j±{x) 
becomes j(x±), where x± are two copies of the spacetime point x: time travels 
from the past to the future for x+ and in the other direction for x_. Sivasubra- 
manian and coll. |44j have proposed to interpret this doubling of spacetime points 
in terms of noncommutative geometry. It would be interesting to follow this track 
for our quadrupling of spacetime points. 



Structure of Green functions 



11 



From the practical point of view, the main applications of our scheme will 
be for the calculation of strongly-correlated systems, in particular for the optical 
response of some materials, such as gemstones, that remain beyond the reach of 
the standard tools of contemporary solid-state physics. 

After the completion of this work, we came across a little known article by 
Sergey Fanchenko, where the cumulants are used to define an effective action [45 . 
His paper is also interesting because it gives a path integral formulation of quantum 
field theory with a general state. His approach and the one of the present paper 
provide complementary tools to attack nonperturbative problems of quantum field 
theory with a general state. 
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